The electron gas with a strong pairing interaction: Three particle correlations and the 
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We derive simplified Faddeev type equations for the three particle T-matrix which are valid in 
the Hubbard model where only electrons with opposite spins interact. Using the approximation 
of dynamical mean field theory these equations are partially solved numerically for the attractive 
Hubbard model. It is shown that the three particle T-matrix contains a term vanishing ~ T^ at 
the Thouless (or BCS) instability where the two-particle T-matrix diverges. Based on the three 
particle term we further derive the low density - strong coupling extension for the two-particle 
vertex function. We therefore understand our equations as a step towards a systematic low density 
expansion from the weak coupling BCS theory towards the stronger coupling limit. 
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I. INTRODUCTION 

One of the key questions which remains in the under- 
standing of the unusual superconductors, like the high- 
Tc-materials, some heavy Fermion compounds or organic 
superconductors is connected with the unusual short co- 
herence length of the Cooper pairs. Even without hav- 
ing clarity on the microscopic pairing mechanism and on 
the evolution of a superconducting phase by e.g. dop- 
ing an antiferromagnetic insulator with holes, this short 
ranged pair formation points towards an effective strong 
electron-electron interaction. It is therefore essential 
to get a deeper understanding especially of the thermo- 
dynamic properties and of finite temperature quantities 
of the electron gas with a strong attractive interaction. 
From the theoretical understanding the electron gas with 
an attractive interaction has - besides the atomic one - 
only two well understood limits. 

The first one is the weak coupling three dimensional 
(3D) case. In this case no two-particle bound state has 
the ability to empty the one-particle continuum and the 
BCS instability can occur for pairs with zero total mo- 
mentum at the Fermi surface. BCS-superconductivity of 
weakly bound pairs with large coherence lengths occurs. 
These pairs simultaneously establish pair formation and 
phase coherence of the BCS wave-function. 

The second limit which is well understood is the two- 
particle problem which gives the zero density limit. Solv- 
ing the two-particle problem in one and two dimensions 
leads for any strength of the attractive interaction to the 
formation of a bound state. Since a pair of electrons can 
be regarded as a boson, the trapping of a pair in the 
bound state can be understood as a "T=0 condensation" 
of a boson into the bound state. Solving the two-particle 
problem in three and more dimensions leads to a clear cri- 
terion to distinguish between the weak and strong cou- 
pling limits. In 3D there exists a critical strength of the 
attractive interaction below which no bound state occurs. 
This marks the weak coupling regime where at finite den- 
sities the BCS instability will occur. Above the critical 
coupling strength the pair of electrons condenses into the 



bound state. 

The simplest model Hamiltonian where the crossover 
from the weak coupling regime towards the stronger cou- 
pling limit can be investigated is the Hubbard model 
with a non retarded, attractive, on-site interaction, even 
though it only shows pairing in the s-wave channel. A 
more complicated correlation including next neighbor 
terms can be included later on where superconducting or- 
der parameters can occur which have non spherical sym- 
metries. Before extending the Hamiltonian we want to 
understand the crossover for a simple on-site correlation. 

An approach has been widely used tQ-.which we refer 
in agreement with many other authorstli3 as T-matrix 
approach. It accounts for both of the limits discussed 
above. It exactly solves the two electron problem in the 
center of mass system and therefore gets exact in any di- 
mension and for any strength of the interaction at zero 
density. The non-self-consistent T-matrix in the 3D weak 
coupling case leads to the same equations to determine 
Tc as the BCS theory. The formation of a BCS wave- 
function can be understood as a macroscopic occupation 
of a bosonic bound state at the Fermi energy by non- 
interacting pairs. The BCS gap-parameter is then pro- 
portional to the density of pairs in the condensate wave 
function. The most problematic part in applying the T- 
matrix to finite densities and to the stronger coupling 
limit is that it neglects pair-pair inteijaiitions completely 
in the non-self-consistent formulationOu of the theory. 

Overall, a self-consistent treatment of the two-particle 
T-matrix seems to compare well with quantum Monte 
Carlo results.B However, certain features, like e.g. a pseu- 
dogap, are either not reporl^ad at all from a self-consistent 
treatment of the T-matrixolH or are not seen in the in- 
tegrated density of states but only in the k-dependent 
spectral functions at extreme parameter ranges.u There- 
fore, it is not tQn_surprising that there is a discussion 
in the literatureEl£3 on which degree of self-consistency 
is the right one, which is also motivated by the correct 
emerging of a superconducting phase at low tempera- 
tures. Using any degree of self-consisitency accounts in 
some way for pair-pair interactions. tlH This is however 



occurring in a nonsystematic way, since a small parame- 
ter is missing. 

In this work we want to tackle the question whether 
there is a systematic way of extending the T-matrix to- 
wards the stronger coupling and finite density regime 
while keeping the complete dynamics of the equations. 
As a first step towards this goal we investigate the three 
particle channel. In the three particle channel the three 
body problem is usually decomposed into two-body frag- 
mentation channels, which define a set of Lippmann- 
Schwinger equations. In general, there are three possi- 
ble two-body fragmentation channels, i.e. three ways of 
combining two out of three particles to interact via a two- 
body force, and a three body interaction. However, for 
the attractive Hubbard model, only two two-body frag- 
mentation channels are left since electrons with identical 
spin do not interact and an ejtplicit three body interac- 
tion is not present. Faddeevtll derived a set of coupled 
equations for the wave functions in momentum space, 
which include the known solution of the two-particle 
problem: the two-particle T-matrix. An equivalent for- 
mulation was obtained for the traasition amplitudes in 
the three body scattering problemJia In this work we fol- 
low the same spirit and calculate the three particle T- 
matrix as an infinite sum over all possible combinations 
of two-particle T-matrices. Since we are interested in a 
dynamic instability, we preserve the complete dynamics 
of the equations m contrast to earlier treatments found 
in the literaturellj'O where one was mainly interested in 
ground state properties. 

In order to be able to do an explicit calculation we 
use k-averaged quantities. Doing so, we arrive at equa- 
tions which on one hand become exact in any dimen- 
sion in the strong coupling limit since in the strong cou- 
pling limit all dispersion disappears. On the other ha 
our equations meet the dynamical mean-field theory,! 
which gets exact at infinite dimensions for any coupling 
strength since the problem can be mapped in infinite di- 
mensions onto a local impurity problem because all cor- 
relations become local. The two-particle T-matrix has 
been studied in the absence of k-dispersionO'El and it has 
been shown that the essential physics of the Thouless in- 
stability remains intact. There is however one important 
difference between the limit of infinite dimensions and the 
strong coupling case. In infinite dimensions one clearly 
has to distinguish between zero momentum and momenta 
which are non zero.Q In the strong coupling limit the k- 
dispersion vanishes. Therefore it is not necessary to treat 
k=0 separately.El 

In the two-body problem the instability in the 
particle-particle channel is called 'Thouless instability'. 
It is given by the equation 



1 + U xiK ^o,n = o) = o, 



(1) 



where U is the attractive on-site interaction and x(-^i ^) 
is the susceptibility in the particle-particle channel. We 
show in this paper that in the three particle channel the 
equivalent of this equation will be given by a determi- 



nant of a matrix which becomes zero at the formation 
of a three particle bound state. It has been discussed 
earlier that three particle terms should be included,tJ 
although they only lead to logarithmic improvement in 
2DE§t3 with respect to a low density expansion. How- 
ever, a systematic low density expansion should contain 
the three particle terms as a next step beyond the two- 
particle T-matrix. In the case of nuclear many body the- 
ory the a-particle as a four particle bound state is of 
partipular low energy and has been discussed by Ropke 
et alM 

The paper is organized as follows: In section || we de- 
rive simplified equations for three particles valid in the 
Hubbard model where electrons with opposite spin do not 
interact. The three parti cles scatter in the vicinity of the 
Fermi surface. In section HI we show haw in a close anal- 
ogy to an earlier work of Gorkov et al.EJ the three particle 
terms can be used to derive a systematic low density - 
strong coupling extension for the two particle vertex of 
a pair of electrons with opposite spins. In section IV we 



apply a k-average and show that the equations reproduce 
the known results in the two-particle channel. We pay 
special attention to keeping the number of Matsubara fre- 
quencies small which are needed for the calculation. The 
analytical extensions to correct the sums over Matsubara 
frequencies are given in appendices. Having established 
a numerical solution of the two-particle problem which 
involves only a minimum of Matsubara frequencies, we 
are able to tackle the three particle instability in section 

0- 



II. THREE PARTICLE INTERACTIONS 



All information on the three particle problem is con- 
tained in a six-time correlation function 



(2) 



Where the cl^,Ci^a denote creation and annihilation op- 
erators of an electron with spin a on the lattice site i at 
imaginary time r, respectively. The brackets (..) denote 
the thermal average and T,- is the time ordering operator. 
In the first glance it seems to be hopeless to calculate such 
a correlation function. However time translational in- 
variance in thermal equilibrium and spatial translational 
invariance on the lattice reduce the number of indepen- 
dent variables to five. Performing all calculations in fre- 
quency and momentum space will further allow via the 
Faddeev equations to solve as well in the "input" as in 
the " output" channel the occurring two-particle problem 
by integrating over the relative momentum and energy of 
a pair. We therefore arrive at functions for the three par- 
ticle channel which depend on three dynamical variables 
only. It is convenient to chose the total energy and mo- 
mentum of the triple as one of them. 

As already mentioned in the introduction, the three 
particle terms can be built up by a repeated use of two- 



particle interactions B3 The special simplification in the 
case of the attractive Hubbard model with only a non 
retarded, on-site attraction, U < 0, is given by the fact 
that first there is no direct three particle interaction and 
second that the correlation U is only effective between 
two charge carriers of different spin which arc on the 
same lattice site. Therefore there is only one nontrivial 
three particle channel where two charge carriers carry a 
spin a and the third one has opposite spin —a. In the 
following, we use a convenient short hand notation for the 
combination of summations over Matsubara frequencies 
and fc-space summations. 



^E/K)-^E^E/(-i'^^) 



(3) 



We use capital Greek letters for Bosonic and small ones 
for Fermionic Matsubara frequencies. Using this nota- 
tion, the Green function of the unperturbed system reads 



G°'^(^i) 



1 



^l-<k^) 



M 



(4) 



where /i is the chemical potential and a refers to the spin 
index which is for our system of electrons (or holes) ±1/2. 
The susceptibility in the particle-particle channel for a 
pair of electrons with opposite spin is 



iyly I 

g A^ N A^ i^i -e(fci) + u 



la'i, 



.i-e(ii:"-fci)+// 



(5) 



The usual two-particle T-matrix results from an infinite 
geometrical series which results from solving the (in this 
case trivially solvable) Bethe-Salpeter equation. 



T{nh) = 



u 



1 + Uxi^h) 



(6) 



This is a vertex function which adds all possible inter- 
actions in the particle-particle channel to the usual di- 
rect Hubbard interaction. It contains the full solution 
of the two-particle problem in the center of mass coor- 
dinate system and therefore allows the formation of a 
two-particle bound state if the denominator 1 + U xi^h) 
becomes zero. This is occurring first at the chemical po- 
tential (iri^ = 0) and for zero total momentum of the 
pair K'^ = and defines the Thouless criterion for the 
superconducting instabilityo Note that we use a differ- 
ent laptation for the sign of xi^h) then in a previous 
worm which is now in agreement with the book of e.g. 
SchriefferEJ but is different from the one used e.g. in Ref. 
Ea. For later use we define the functions 






(7) 



with I S {0, 1, 2}. Here r*^°)(r2;j) is a two-particle propa- 
gator and r'^'^\tlh) is the above defined vertex (g) with- 
out the single interaction between the particles. The def- 
inition of the different T^^^Qh) will be needed later on 
to be able e.g. to treat the contribution of the Hartree 
term separately. ■— ■ 

Following the ideas of Faddeevllil, the three particle 
terms can be broken up into a two-particle problem plus 
an additional particle. Therefore, we define the following 
functions. 

Fi(wi,c^2,W3'):=-G°"(c^i)G°""(a;2)r(c^2+c^3') (8) 

This function is illustrated in Fig. |l|a. It contains all 
possible repeated correlations between particle (2) and 
(3') and contains the propagator for particle (1) which 
does not interact with any of the other two particles. F 
is already a part of a fully interacting six-time correlation 
function. However the time and spatial translational in- 
variance together with solving the two particle problem 
in it's center of mass system has in frequency and momen- 
tum space reduced the function to an expression which 
contains only three independent dynamical variables. 
The second function which we define is: 

F3{uj3,UJ2',UJi) ■= -G^^iuiy) G° ''{UJ2') r(LJi +UJ2') 

(9) 

This function is shown in Fig. Hp. The functions -Fi 
and F3 represent the only possible two-body fragmenta- 
tion channels which are present in the attractive Hub- 
bard model. A repeated application of F^ and Fi turns 
the frequency UJ2> into an internal frequency over which 
a summation has to be performed. A diagrammatic il- 
lustration of this term is shown in Fig. 13. 

-^F3{uj3,uj2',uji) Fi{uji,u;2,(^3') S^2+^y,uj2'+^3 (10) 

^ 2' 

The delta function stems from energy and momentum 
conservation for a pair. It is more convenient to intro- 
duce a new variable. This is the total energy (and mo- 
mentum) of all three particles cug = uji+uj2+ ^3' since it 
is a conserved quantity for the three particle scattering 
events. In this way we define 

M^^'>{ujg,uj3,0J3,) := 

■^ E ^3 (^3' , W2' , Ci^g -UJ2' -1^3') 
^ 2' 

X Fi{ujg - LU2' - W3',LJ2' + WS' -^3,^3). 



fll) 



Note that all the Pauli blocking factors which are of im- 
portance for t he . di scussion of three body interactions in 
nuclear mattercSEZI are included in the summation over 
Matsubara frequencies due to the analytic properties of 
thermal Green functions. 



We can now successively construct all possible interac- 
tions among three particles by applying alternate scatter- 
ing in the only possible two-body fragmentation channels 
by a repeated application of Ai^^' . For example applying 
M^^' twice gives one additional internal frequency (and 
momentum) and therefore yields an additional summa- 
tion. 

Af (2) (cc^, 0,3,^3,):= 

' 3" 

= :M(^)®m(^) (12) 

This defines a formal product for the alternate interac- 
tion. In appendix |A| we show that this formal product 
can also be written by introducing a metric. This prod- 
uct can be repeated several times: 



M(")(L^g,C^3,C^3') 



1 



/?(" 



X! M'^^Hujg,UJ3,UJ3n 



3",3"',...,3(") 

X M(^)(wg,u;3",W3"') ... M^^^ujg,uj^(„),ujy) 
= Af(^) ®M(^) ® .... ® m(^) (13) 

In the same spirit one can define a neutral element for 
the product defined in Eq. (n2h. 



m(") (E)E = m(") 



(14) 



which leads to: 



It can be seen immediately that the function E fulfills 
all properties of a neutral element with respect to the 
product defined in Eq. (12). 



M^^'^^E = i^M(")(c.g,c.3,c.3")/3<5.,3,„..3, 



M^"'>{ujg,UJ3,UJ3') 



(16) 



We are now able to express all possible interactions which 
can occur between three particles in the attractive Hub- 
bard model. This is given by the infinite sum over all the 
three body interaction terms M^"\ It yields a geometri- 
cal series which has to be calculated. 



which can be formally written as: 

M'°*(u;g,u;3,c^3')= \p S^,,^^, ^ M^'\u;g,u;3,t03') 



The inversion has to be such that 



M'°*(g) 



E - M(i) 



E. 



(18) 
(19) 



If on decreasing temperature the function E—M'^'^^ ceases 
to be invcrtiblc this defines an instability in the three par- 
ticle channel. Such an instability marks the appearance 
of a three particle bound state. In section ^ we show that 
in dynamical mean-field theory, which means in the ab- 
sence of fc-dependence where iOh — > ia;„, the condition 
for the instability can be expressed as 



det[£:k-. 



Af, 



(1) 



J-0- 



(20) 



This determinant has to be calculated for an infinite num- 
ber of infinite dimensional arrays. However, restricting 
the number of Matsubara frequencies to a finite number 
allows to calculate the condition for an instability in the 
three particle channel. This is the three particle equiva- 
lent of the Thouless instability (o) which occurs for scat- 
tering of two particles. We like to mention that a general 
formulation for an arbitrary interaction has been given 
by Ethofcr and Schucko Their equations in the p-p-p- 
channcl contain the complete dynamics as our present 
work does. Due to the particular simple interaction of 
the Hubbard model, our equations obtain the much sim- 
pler form given in thi'A article. 

In nuclear physicsE£IEZI it is often convenient to use sim- 
plified formulas for the three-particle T-matrix which de- 
pend only on one dynamical variable (z) instead of an 
expression which depends on three dynamical variables 
{ujg,uj3,oj3i). The reason for this is that usually only 



the dependence on the center of mass {z 



j) is con- 



sidered. Such equations are recovered when all internal 
initial times are set to zero when therefore all internal 
dynamics of the triple is neglected: 

A/(wg,k3,k3/) = liin iuj^ iuj^, M{ujg,uj3,LUy) (21) 



In section IV 



we restrict ourself to an approximation 
which, as we will show, preserves the main physics and 
mathematics of the two particle Thouless instability to- 
gether with the complete dynamics of the equations. Af- 
terwards, we are going to show how the physics is altered 
due to an inclusion of the three particle terms. 



M'°\uJg,CU3,UJ3') = ^A/(")(c^g,C^3,^3') 



n=0 



fiS, 



UlS.UJ^I 



M^'^\ujg,U3,uj3>) 



-^M(l)(Wg,W3,W3")A^''H^ff,^3",W3') 
P 3" 



(17) 



III. SYSTEMATIC LOW DENSITY - STRONG 

COUPLING EXTENSION OF THE TWO 

PARTICLE T-MATRIX 

In this section we show how the three particle T-matrix 
can be incorporated to derive a low density strong cou- 
pling extension of the equations in the particle-particle 



channel which lead to the BCS instability. For a contin- 
uous system with a non retarded interaction it has been 
shown by Gorkov and Melik-BarkhudaevEJ that a term 
which is first order in low density [^ kpa) can be used 
to obtain a low density correction for the two particle 
instability. The equations ofcll contain however only a 
correction which is of second order perturbation theooz. 
in the interaction. The term which was considered intiJ 
is illustrated in Fig. ^. In the Hubbard model this cor- 
rection is exactly zero since a direct interaction between 
particles with the same spin is not present. To calculate 
the next low density correction for the attractive Hub- 
bard model the first non trivial term is of third order in 
the interaction U (the term of second order will only lead 
to a Hartree shift). It is illustrated in Fig. ^ 

U°ein,^i,T)xm^Uxm + 



into the Thouless criterion: 



i + [/eff(J7 = o,r,/i)x(f^ = o) = o 



(26) 



The correction due to the three particle terms vanishes in 
the zero density limit. In the limit of vanishing density 
the summation over l/fij^s' [H (E) AI*°^ gives zero and 
Eq. (J^^ gets identical with the instability of the two 
particle T-matrix. The correction due to the product 
[H (g) M*°*] vanishes also at high temperatures since all 
pair propagators go to zero at infinite temperature 



lim Y(f7) = 



(27) 



/33 



2,3,4 



G^{u^)G^{^ + Ul + LOi - UJ2 - LOs) 



(22) 



The infinite sum over all the terms up to infinite order 
of the interaction U involves the three particle T-matrix 
M*°*(u;g,a;3, wa/) and is shown in Fig. |^. To obtain this 
we define the function: 

iJ(f) - 6^3', ^3', C^3) = ^ E ^"(^2)G"(f} - C^a) 
' 2 

r(c^3. + c^2)G°(c^2 + c^3' - ^3)G"(a;3) (23) 

which can be obtained by replacing one T-matrix in the 
derivation of Eq. (O) by the bare interaction U. There- 
fore one can obtain the effective two particle interaction 
by calculating the product H (g) M*°*: 



Eq. ( pq ) gets further exact in the weak coupling regime 
where it reduces for all densities to the Thouless or BCS 
condition for a two particle instability since the lowest 
order correction to U is of third order in perturbation 
theory. 

A numerical solution of Eq. (24) for the limit of van- 
ishing fc-dispersion is in progress and will be published 
elsewhereEa. It is also possible to solve Eq. (g^ in a self- 
consistent way where selfconsistency has to be achieved 
for the function t/cff(17, u, T). This might enlarge the 
applicability of Eq. (gj) even beyond the low density 
regime. We further like to mention that the freque ncy 
and momentum dependence of C/cff(f^, /^,T) in Eq. (|24| ) 
does not lead to a retarded interaction since the variable 
ri refers to the total momentum and total energy of a 
pair who's interaction is being treated. 



IV. SOLUTION OF THE TWO-PARTICLE 
PROBLEM 



u,sin,^,,T)xm^Uxm 



^ 3' 

2^/2 



/9 



We need a reliable solution of the two particle prob- 
lem, if we want to solve the three-particle problem since 
the two-particle T-matrix enters into the Faddeev equa- 
tions. We will solve the three particle problem only in a 
k-averaged approximation, therefore we show in this sec- 
-, tion that (i) we are able to reproduce the known results 

2_^G^{Cl — 012)0" (0^2) G° (0^2) -r yjG"(i2|) of the two-particle T-matrix using a k-average and (ii) 

we demonstrate that we are able to solve this equations 
using a minimum of Matsubara frequencies which is the 
essential basis to calculate in Sec. ^ quantities for the 
three particle problem. 

In the noninteracting limit U — the free k-averaged 
Green function is 



(3 



where the last term is the term of second order which 
contains the Hartree term. It has to be subtracted in or- 
der to obtain a consistent inclusion of the Hartree term. 
The summation over ws' closes all the possible three par- 
ticle terms with a hole line making Eq. (p3) a low density 
strong coupling correction. The factor of 2 stems from 
the two possible realizations of the three particle correc- 
tion which are tjt and |ti ^s combinations for the spins 
of the triple. The product is given as: 

[7?®M*°*](Wg,W3,^3') 

= -^i/(a;g,W3,W4)M*°*(wg,c^4,c^3')(25) 

l-^^ 4 

Equation (g^ leads to the low density - strong coupling 
extension of the two particle instability when inserted 



G"(iuj„ 



Die) 



_^ lUJn -CJ + fi 



-de. 



(28) 



D is the noninteracting global density of states. Since we 
are heading towards a systematic low-density expansion 
we are mainly interested in low band fillings. Therefore 
we approximate the density of states I? by a rectangular 
shaped one, which shows the correct analytical behavior 
for a 2D system at the lower and upper band edge. 

Die)^l^[Qie + W/2)-eie-W/2)] (29) 



W is the bandwidth. This density of states results in 



1 



G'"(iw„) = —[HiLUn + W/2 + ^) 



W 



lll{iUn-W/2 + ^J,)]. 



(30) 



The k-averaged pair susceptibihty can be calculated from 
the k-averaged Green functions. 



K 



(3N^ 



Y^ G°ik,iLu,n)G"iK ^k,in,, 



K .k.m 



i^GO(*c^„0GO(*f]„-ic^™) 



(31) 



The anal ytica l correction for calculating x is shown in 
appendix B 2 . We approximate the fc-averaged vertex 
function r(if2„) neglecting mean squared fluctuations of 
the pair susceptibility.^ 



r^^H^^n) 



1 + Ux{l^n) 



(32) 



r^^) {i^n) is used (see Eq. (0)) which contains the infinite 
sum of all possible interactions within the pair except the 
single interaction. The single interaction U only leads to 
a Hartree term when inserted into the selfenergy The 
Hartree term shall be considered separately. The full 
interacting Green function is 



with the self-energy 



(33) 



l]0(zw„)-- Y. r(2)(jc^„-fic^„)G°(zc^„,). (34) 



1. Particle number n 

Since we need the two particle terms in the three par- 
ticle T-matrix we demonstrate in the following that we 
are able to obtain the main results of the two particle 
T-matrix calculation by applying the k-average. We put 
special emphasis on keeping the number of Matsubara 
frequencies small. 

The expectation value of the particle number with a 
certain spin a is calculated from the single particle Green 
function. 

(n.)(/3,/i) = lim- Y G{^coly'^"' 



-^ReG(zc.,) 
'^^ 1=0 



(35) 



The total particle number n is the sum over the two spins. 



(n) = 2{n,) 
To obtain Fig. W^ we have solved Eqs. (t50| 



(36) 

_ (H non- 

self-consistently. We have plotted the lines of constant 
densities (n-f ) — (ni) — {n)/2 in the plane of the thermo- 
dynamic parameters, the chemical potential without the 
Hartree contribution jl = jj, — U{n)/2 and the tempera- 
ture T. We were able to use only a maximum number of 
-^max = 60 Matsubara frequencies by applying the ana- 
lytical corrections which are shown in appendix |^. Using 
these Nyaax Matsubara frequencies we reach temperatures 
which are as low as ksT — 0.01 [W/2]. One has to assure 
that the range which is covered by the A^max Matsubara 
frequencies is larger than the bandwidth of the system 



,(2iV„ 



l)7r 



/3 



>W 



(37) 



In this way less than one hour computer time on a Pen- 



tium II with 450 MHz was needed to obtain figure [V] 



All lines of constant density collapse into the bound state 
level on lowering the temperature. A result which has 
been first obtained by Schmitt-Rink et al.u who used a 
non conserving approach. In our approach the particle 
number is conserved and reaches 1 at the Thouless insta- 
bility. 



2. Double occupancy 

In this section the two particle correlation {n-^ni^) is 
computed from the two particle propagator F'"' . 



Wni) 



1 



E 



r^°\^n^) 



(38) 



In Fig. VI we have plotted lines of constant double occu- 



pancy (n^ni) in the plane of the thermodynamical vari- 
ables /i,T. We have chosen the values of (n-fni) such 
that some match the square of {n)/2 in Fig. N^. In 
the limit of high temperatures or of vanishing interac- 
tion (U — > 0) all correlations disappear and (ni^ni) = 
((n)/2)^ To be specific: At high temperatures the line 
with (ntni) = 0.04 will join the line with (n) = 0.4 in 
Fig. ^ and the one with (n-fn^) = 0.09 will join the 
line with (n) = 0.6. On lowering the temperature how- 
ever, correlations become important and (n]'ni) exceeds 
the value of ((n)/2)^. Close to the Thouless instabil- 
ity (n|n|) even diverges. A diverging value of (rifn^) 
marks an unphysical behavior since a pair of particles 
cannot be more then completely correlated. Therefore 
{n^ni) should not exceed the value of {n)/2. The behav- 
ior in the non-self-consistent T-matrix calculation can be 
understood since summing over bosonic Matsubara fre- 
quencies in Eq. ( ^ occupies the bound state with a 
Bose distribution and (7i|?i|) — > 00 marks the macro- 
scopic occupation of the bosonic bound state with non- 
interacting pairs. A large part of this unphysical behav- 
ior is due to the neglect of pair-pair interactions in the 



non-self-consistent T-matrix. Self-consistency seems to 
repair such a behavior but is uncontrolled since a small 
parameter is missing. In order to achieve a systematic 
expansion with finally the density as a small parameter 
we perform a numerical investigation of the three particle 
terms in the next section to answer the question, whether 
this unphysical divergence could be canceled by the sum 
of all three particle scattering processes which contribute 
to {n-^ni). 



V. THREE PARTICLE CHANNEL 



In the k-averaged approximation the three particle 



renormalization is caused by a third particle as derived 
in Eq. (||). 



terms are obtained by replacing ujk 
we obtain 



IUJ„ 



M'°\iujllCOrn,iuJm') = \ E - M '■'\lU?,) 



Doing so, 



(39) 



for the three particle ladder M*°'. This means that in 
the k-averaged approximation M*°' can be obtained by 
inverting the right hand side of Eq. (B9h and multiplying 
it with a prefactor of /3^. However, this is an infinite (due 
to iuj^) number of inversions of infinite dimensional ma- 
trices which has to be performed. Our efforts described 
in the previous section |^ and in appendix ^ enable us to 
reduce the number of Matsubara frequencies needed for 
such an inversion to such a low value that the inversion 
of Eq. (p9) can be performed. In the present work we 
want to answer the question what is the behavior of the 
three particle terms at the two-particle instability. 

In Fig. H we have plotted the logarithm of the deter- 
minant 



det pS, 



M^^\i 



■A.'^ 



(40) 



for one chosen total Matsubara frequency icj^ — it: / (3^ 
for an attractive interaction U ~ —2\W/2\ at a chemi- 
cal potential jl — — 0.5[W^/2] for which Tc turns out to 
be Tc « 0.33301 [W^/2]. On lowering the temperature the 
determinant diverges at the Thouless instability propor- 
tional to (T - Tc)-2. This is plotted in Fig. | for 60 and 
120 as the maximum number of Matsubara frequencies. 
In contrast to this, the equivalent of the determinant in 
the particle-particle channel (see Eq. (|l|)) goes to zero 
linearly in (T — T^) at the instability. This has two im- 
portant consequences. The first one is that even with 
a direct inclusion of the three particle terms the double 
occupancy will still diverge at the Thouless instability, 
this is due to the fact that the double occupancy will be 
calculated as a sum of all two particle terms which di- 
verge at the Thouless instability plus contributions from 
the three particle terms which vanish at the two-particle 
instability and have therefore no chance to cancel the 
diverging two-particle terms. The second consequence 
is therefore that a correct inclusion of the three particle 
terms has to go via a renormalization (possibly even in 
a selfconsistent way) of the two particle instability. The 



VI. CONCLUSION 

In this article we present simplified equations for the 
three particle terms of the attractive Hubbard model. 
In particular the fact that in the Hubbard model elec- 
trons with identical spin do not interact directly with 
each other, enables us to simplify the equations consid- 
erably^ Based on the three particle terms we derive in 
sec. [II a low density strong coupling extension of the 



two particle vertex function (Eq. (24)). A third particle 
renormalizes the effective interaction in the two particle 
channel and therefore leads to a density and coupling 
strength dependent reduction of Tc- 

In order to solve some of the equations we derived, 
we applied a k-averaged approximation which becomes 
formally equivalent to the dynamical mean-field theory 
where all correlations become local. Using such a k- 
averaged approximation we are able to actually calcu- 
late three particle quantities with its complete dynam- 
ics. We investigate the three particle ladder, where the 
condition for the Thouless instability in the two-particle 
channel (sec Eq. (1^)) is now replaced by a determinant 
of a matrix. This determinant diverges proportional to 
(T — Tc)^^ as the Thouless instability is reached. In the 
near future it should be possible to extend the present 
calculations in three ways. First, the double occupancy 
as a function of the thermodynamic parameters includ- 
ing all three particle terms should be investigated. Cal- 
culating this quantity at low densities should answer the 
question which degree of self-consistency is the best one if 
the two-particle T-matrix is evaluated self-consistentbi.. 
This question is currently discussed in the liter ature.EIEHI 
Second, at temperatures below the Thouless instability 
one should continue the determinant (see Eq. (^)) to 
the real axis with respect to iw^ in order to obtain a 
condition for a three particle bound state which should 
arise as the determinant vanishes at the real axis. The 
most important extension is a numerical exploration of 
Eq. (^^ in the space of the thermodynamic variables /i 
and T. This is currently under investigation. 
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APPENDIX A: PRODUCT OF THREE PARTICLE 
INTERACTIONS INVOLVING A METRIC 



The formal product defined in Eq. (|fj) can also be 
written by introducing a metric tensor 



g- 



(3 ^v = P S^^ .^iN6^,M, 



(Af) 



where the last term reminds on the definition in Eq. (0). 
This leads to a co- and contravariant notation for the 
matrices using the Einstein sum convention. 



M: 



33' 



The formal product then reads: 



M(g)N = M33,, TV^ 3, = M33,, 5-3"3"' ^^,„^, 



(A2) 



(A3) 



fn the derivation presented in se c. W]\ and [II we did only 
use the notation defined in Eq. ([l^ 



APPENDIX B: ANALYTICAL CORRECTIONS 

The infinite sums over Matsubara frequencies cannot 
be evaluated nuinerically. iVmax denotes the maximum 
number of Matsubara frequencies that we want to sum 
over numerically. To preserve the features of the infinite 
sums, the addends are approximated for large Matsubara 
frequencies by terms, which can be summed analytically 
for an infinite number of Matsubara frequencies. 



1. The particle number 

The real parts of the interacting Green functions 
in (Bq) are approximated by the noninteracting Green 
functions for large Matsubara frequencies (MFs). The 
noninteracting Green functions exhibit a l/(iw„ + /i)- 
dependence for large MFs. 



G°{iLij„) ""-14"'''° Gs{itOn) 



1 



(Bl) 



have to be approximated by Gg for values of iw„ which 
are far away from the real axis, e.g. for n > iVmax Sind 
n < -iV,n 



1. 



XS 



(i^m) = ^ X! Gs{iuJn)Gs{i^m- i(^n) 



tanh^ 



(B3) 



2rmr {2n + l)Tr 
[2(m-n) - l]7r 

^(m-n-l) 



2N„ 



Xii^r. 



-g ^ [Gh{iu]n)Gh{ii^(m-n-l) 

^ ra=2Ar„ax-l 

-GsiioJn)Gsii^m - Jw„)] 

+Xsii^m) 



(B4) 



with 

GhiiuJi) 



G°iiLUl), -^,„ax - 1 < / < iV„ 

Gs{iuJi), otherwise 



(B5) 



In subsequent chapters it will be necessary to ap- 
proximate the pair susceptibility for large Bosonic MFs. 
Therefore we replace the noninteracting Green functions 
in (|3l|) by their definitions (p8|). This results in two inte- 
grations over e and e'. 



X(«^m) 



1 



f3W^ 



E 



n— — 00 



de 



de'- 



1 



i{flrn - tJn) - c' + /i 



iuJn - e + l-J. 

(B6) 



N„ 



n{fi,ij) 



^ Y, ReG(zu;,) -f | ^ ReGs{iiOi) + \ 



-y 



a 



i = iV„,ax + l 



ReG(ia;0 - -^ 



A^/5 , 1 



(B2) 



The 1/2 stems from the sum ovetthe imaginary part 
where a sum over sin(a;)/x occurs. Ej 



2. The pair susceptibility 

For the numerical calculation of the pair susceptibility, 
we look at the infinite sum (^ll) . The two Green functions 



If we interchange the two integrations with the sum- 
mation over TO, carry out the summation, transform the 
integration coordinates to x = e -I- e' and j/ = e — e', we 
can integrate over y analytically. The resulting spectral 
representation of % can be approximated for large Mat- 
subara frequencies by xw ■ 



/w(/3,/i) 

Xwii^ra) 



1 ^^ l + cosh[/3(I¥/2-^)] ^g^^ 



m 1 + cosh[/5(VF/2 + ^)] 

iflm + 2/i 



(B8) 



Here /^y denotes the averaged spectral function for x- 
The function xw has been used to obtain the results 



shown in section IV 1 and IV 2 



3. The self-energy 



to correct the sum in Eq. (t58) analytically. 



Looking at the definition of the self-energy (M) we 
have to approximate the vertex F*-^'. It can be approxi- 
mated in first order of j^ for large Vtm by U'^xw- Hence 
we calculate the following infinite sum. 



-J CXj 



a 



ibJm + M 



o/3m 



-/3m 



(B9) 



The analytically corrected sum over a finite number of 
MFs is a term of the form 



S(iWm) 



2N„ 



y^ G° {iuJn)'^^'^\ii^m + «W„) 






(BIO) 



It gives the expression of Y,{iujm) which contains the nu- 
merically summed expression for the selfenergy for Mat- 
subara frequencies close to the real axis and contains the 
analytical correction for MFs far away from the real axis. 



4. The double occupancy 

The two particle propagator T'^^\i^m) equals 
Xw{i^m) at large f2m- Therefore we utilize 



771 — — C30 



Xw(«^m) 



fwiP,fJ-) 

I _ g-2/3M 



(Bll) 



5. Three particle ladder segment M^^' 

For the analytical correction of the sum in Eq. dlQ ) 
again for the case of vanishing k-dispersion we look at 
the UJ2' depending terms. Here the terms of the product 
in Eq. (uw which do not contain uj2' are not shown. 



i X) Gsiiuj2' +iuJ3-iuj3')Gs{i(^g 



IUJ2' 



IUJ3) 



xGs{iuj2') [~U + U\w{i^3 + i^2')] (B12) 
-U [d%bb' + U {bb' + aoo + a^bb') fw{P, ^i)\ 



with 



aabbb'cc' 



Go 

a 
a 
b 
b 

b' 
c 

c' 
d 
o 



- ^/3m 



= e'^'^ + l 

= e'"^-! 

— 3^ + iujg 

= ^ + iuz 

~ fi + iuj3> 

= 2/i — iuj3 + iLUg 

= 2/i — iuj3' + iLOg 

= iuj3 + iuj^i 

= w^ -f IJJ3UJ31 — iujg(d — 4/i) — Adfi + 2/i^ 



(B13) 



This allows to calculate the three particle function 
M{Ljg, ijj3, uj^t) for the case of a neglected k-dispersion. In 
section ^ the determinant of E — M close to the Thouless 
instability is investigated. 
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. tJl, 



■'^2, 



•^l, 



. CO')/, — O" 



■ ^3, 



FIG. 1: The function Fi(loi,uj2,uj3') is shown diagrammati- 
cally in Fig. hla . The shaded area denotes the T-matrix and 
represents all the possible interactions between particle 2 and 
3. Fig. hb shows the function F3{(jJ:i,ijj2',LUi). 



12 



W, =U)i+U)2+ ^^3' 



wi, a 



W2, -O" 



W3, (T 



UJ3,, a 



FIG. 2: The repeated application of Fi and F3 gives the fiinc- 
tion M^^' {ujg,u>3,Lij3i). This is shown as a diagrammatic rep- 
resentation. 
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FIG. 3: Low density correction which was considered incil. 
The last term does not give any contribution for the attractive 
Hubbard model since a direct interaction between particles of 
identical spin is not present. 
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FIG. 4: First non-vanishing term for a systematic low density 
correction which is of third order in the interaction U. 
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^«=n 





FIG. 5: Infinite sum over all orders of the interaction U which 
gives the complete low-density strong coupling extension for- 
mulated in Eq. (M) with the exception of the Hartree term 
which has to be subtracted. 
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FIG. 6: Lines of fixed particle numbers on the T-p-plane. 
The attractive interaction U is equal to the bare bandwidth 
W and n — 1 denote half filling. At low temperatures all 
lines of constant density collapse into the two-particle bound 
state. 
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FIG. 7: Lines of fixed (n-^ni) on the T-/i-plane. At the Thou- 
less instability (n|n|) divergences. The attractive interaction 
U is equal to the bare bandwidth W . 
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FIG. 8: The determinant given in eq. (U(j) diverges at the 
Thouless instability proprtional to {T — Tc)^'^. This is shown 
for jl = —0.5 and Tc ~ 0.37 for two maximum numbers of 
Matsubara frequencies M — 60 and M = 120. The total fre- 
quency iujn was choosen at n = 1. The attractive interaction 
U is equal to the bare bandwidth W. For a comparison the 
value of 1 + U X which goes to zero linearly in T, with {T — Tc), 
is plotted. 



